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Abstract 

We investigate some properties of a system of Dirac fermions in 



Qh! 2 + 1 dimensions, with a space dependent mass having domain wall 

I like defects. These defects are defined by the loci of the points where 

the mass changes sign. In general, they will be curves lying on the 
spatial plane. We show how to treat the dynamics of the fermions 
in such a way that the existence of localized fermionic zero modes 
on the defects is transparent. Moreover, effects due to the higher, 
no-zero modes, can be quantitatively studied. We also consider the 
relevance of the profile of the mass near the region where it changes 
sign. Finally, we apply our general results to the calculation of the 
induced fermionic current, in the linear response approximation, in 
the presence of an external electric field and defects. 
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1 Introduction 



In this paper we study some properties of a system of Dirac fermions in 2 + 1 
dimensions with a space dependent mass that corresponds to domain wall 
like defects. Dirac fermions models in 2 + 1 dimensions appear as effective 
theories for various systems of interest in particular in condensed matter 
physics. Consider for instance a tight binding model for spinless electrons on 
a two dimensional square lattice, in the presence of an external magnetic field 
of half of a magnetic flux quanta per plaquette. The nearest neighbor hopping 
amplitude is one, the second neighbor hopping amplitude is t/4, and there is 
a staggered potential fi(—l) Xl+X2 , where X\ and x 2 are the coordinates on the 
plane. It has been shown 0] that, in the continuum limit, the effective theory 
for this system corresponds to two massive Dirac fermions whose masses are 
proportional to (/i ± t). In the low energy sector, the heaviest particle can 
be neglected and the theory to be considered is a Dirac fermion with mass 
proportional to M — \i — t. Now imagine that in the system there is a 
domain wall (or stacking fault) created in the process of crystal growth by 
some specially prepared circumstances. At the level of the effective Dirac 
theory this can be obtained by replacing M by — M in half of the plane 
Thus the mass becomes position dependent with a profile M(xi,x 2 ) = 
M(9(x 2 ) — 9(—x 2 )) for a sharp wall at x 2 = 0. Exactly the same effect can 
be obtained when the system, even without any defect, has an odd number 
of lattice rows in one direction. If the system is build up periodically the 
presence of the extra row produces the same effect as a domain wall. 

In this work we study the general case in which the mass of the Dirac 
fermions changes sign along an arbitrary curve in the spatial plane, and with 
an arbitrary profile near the regions where the sign change is observed. Due 
to the Callan and Harvey mechanism [|J, localized chiral zero modes appear 
on the defects. This property has been extensively applied, following an idea 
by Kaplan |J , to the problem of defining a lattice action for chiral fermions, 
thus providing a possible way of overcoming the kinematical obstruction 
posed by the Nielsen- Ninomiya theorem ||. The proposal evolved into the 
somewhat more abstract, but nevertheless equivalent, 'overlap formalism' J7|, 
0. |9|. In Kaplan's formulation, one deals with step-like rectilinear defects for 
the mass in an odd number of spacetime dimensions. Based on those ideas 
and techniques, we shall here study the contribution to the effective action 
of the zero modes appearing along a curve of arbitrary shape in the spatial 
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plane in 2 + 1 dimensions, quantifying the effect of the extended bulk states 
and of the particular shape of the domain wall profile on the dynamics of 
the modes. Localized non chiral zero modes are also taken into account. As 
an example of an application, we use the linear response approximation to 
derive an expression for the induced current in the presence of an external 
uniform electric field. 

The paper is organized as follows: In section 2 we introduce an eigenmode 
expansion appropriate to the presence of domain wall defects, considering 
first the simplest case of a rectilinear domain wall, and then extending it to 
the more general situation of a defect of arbitrary shape. We apply this to 
the particular example of a circular defect. In section 3, the eigenmodes ex- 
pansion is applied to the problem of evaluating the vacuum fermionic current 
in the presence of an external field, to understand the consequences of having 
defects to the spatial distribution of the fermionic current. The results for 
the concrete example of the induced current due to a constant electric field 
are presented and discussed. Section 4 is devoted to our conclusions. 

2 Eigenmodes expansions in the presence of 
domain wall defects 

A functional integral over fermionic (Grassmann) fields usually carries a 
quadratic action, in terms of some operator, which depends on the exter- 
nal fields. They are regarded as sources at this stage, either because they 
have no dynamics (as in our case), or because it would be desirable to post- 
pone their quantisation. The formal integration over spacetime dependent 
Grassmann fields may be converted into one over the (Grassmannian) coef- 
ficients for the expansion of the fields in a complete set of eigenfunctions of 
some operator. A judicious choice of eigenfunctions may simplify the task of 
evaluating the functional integral, or provide a better starting point for an 
approximate treatment. 

In the present case, we deal with Dirac fermions in 2 + 1 Euclidean di- 
mensions, in the presence of an external Abelian gauge field A^, and a space 
dependent parity breaking mass M(x), with x = (xi,x 2 ). The Euclidean 
action S is defined by 

S = fd 3 x^(x) [fl + ieJL + M(x)] (1) 
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where x = (x, £3), X3 being the Euclidean time coordinate. Hermitian 7- 
matrices are in an irreducible (2 x 2) representation of the Euclidean Dirac 
algebra: 

{ 7m ,7,} =2^. (2) 

By definition, the mass M(x) has a domain wall defect whenever it 
changes sign. As these defects proceed from the regions where the surface 
spanned by the function M(x) crosses the plane M(x) = 0, we see that they 
shall consist of (possibly disconnected) curves. Let 

r -> x(t) , r e [ti,t 2 ] (3) 

be a parametrisation of one of the connected components of the curve. Then, 
of course we shall have 

M{x{t)) = 0,Vr E [n,r 2 ] . (4) 

This curve can in principle have an arbitrarily complicated structure, depend- 
ing of course on the assumptions made about the function M(x). We shall 
assume that the mass is regular enough as to define a curve which is a dif- 
ferentiable manifold. We will see that the phenomenon of localization in the 
general case may be understood in a simple way, if one defines a convenient 
system of coordinates in the neighborhood of each point on the curve. For 
each coordinate patch, the physics can be shown to be equivalent to the one 
of a mass term depending on only one coordinate. Thus, in each coordinate 
patch, the defect will look like rectilinear. Our strategy will be to consider 
first the simpler situation of a rectilinear defect, and then to generalize it to 
the case of a wall defining an arbitrary curve, through the use of the above 
mentioned system of coordinates. 



2.1 Rectilinear defect 

This is indeed the simplest situation, consisting of the mass being a function 
of only one of the spatial coordinates, say X2 [], and changing sign along the 
straight line 22 = 0: 

S = J d 3 x ^(x)V^(x) (5) 

1 For a step-like defect, this is the 2 + 1 analog of the configurations relevant to Kaplan's 
formulation Q. 
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where we defined 

V = 7 M D„ + M{x 2 ) . (6) 

and D^ = d^ + ieA^. 

It is easy to show that in order to be able to disentangle the dynamics of 
the fermions into two pieces: one depending on x 2 only, and the other one 
with support in the plane x 2 = 0, the most general gauge field configuration 
we can consider is 

F 2M = 0,^ = 1,3, (7) 

which corresponds to having no magnetic field (F 2 i = 0), and no electric 
field perpendicular to the defect (F23 = E 2 = 0). These are gauge- invariant 
conditions, which in the gauge A 2 = (which we shall adopt) leads to 

d 2 A^ = , n = 1, 3 . (8) 

In this gauge, we immediately see that the operator T> may be rewritten as 
follows: 

V = (a+ f)V L + (a f + i)V R (9) 
where a, a) are operators acting on functions of x 2 , 

a = d 2 + M(x 2 ) , a f = -<9 2 + M(x 2 ), (10) 

Vl,r are projectors along the eigenspaces of the matrix 7 2 : 

V L = ^(1 + 72) , V R = i(l- 7 2), (11) 

and fi is the two dimensional Euclidean Dirac operator corresponding to the 
two coordinates x\ and X3, which we denote collectively by x, namely 

i = li{di + ieA l {x)) + 73 (<9 3 + ieA 3 (x)) . (12) 

Expression @ suggests the possibility of getting rid of the dependence 
on x 2 for the fields, by a suitable expansion in the modes of some operator, 
and obtaining in a way a 'dimensional reduction' from the three dimensional 
spacetime to the two dimensional one corresponding to x. As T> itself is not 
Hermitian, we use instead the positive Hermitian operator "H 

n = v ] v , (13) 
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which, by using the explicit expression (|9|) for V leads to 



H = {h- f)V L + (h- f)V R , (14) 

with h = a) a and h = aa) . The appearance of the two conjugate Hermitian 
operators h, h suggests an expansion for the x 2 dependence of the fermionic 
fields in terms of their eigenfunctions, namely 

ty(x 2 ,x) = [fini^tpP (£) + <f>n(x 2 )lpR ] (x) 

n 

*(x 2} x) = E[* ) (^)4(^)+A n) (£)4(^)] (15) 

n 

where the subscripts L, R denote the 'chirality' defined by the matrix 72 and 

h(p n (x 2 ) = \ 2 n <fr n (x2) , H n (x 2 ) = A^0 n (x 2 ) 

(4>n\4>m) = 5n,m , (4>n\4>m) = $n,m ■ (16) 

To write eq. (|16|) we have used that ^^(5) = VL,R^ n \ and ^%,(x] 



J L,R\ 

iP^Vr^l- Note that the two dimensional fermionic fields ip^ and their Dirac 
adjoints are independent variables for different values of the discrete index 
n. We have made explicit the property that h and h are positive (the X n are 
assumed to be real and we fix their sign, by convention, to be positive) and 
have the same spectrum, with the only possible exception of the zero modes 
A n = 0, since, for any (f> n with X n ^ 0, there also exists one eigenvector of h 
with identical eigenvalue 

Hn{x 2 ) = X 2 n <p n {x 2 ) Jl [rr-a^nM] = [— a<f) n (x 2 )) (17) 

where the factor ^- is introduced to normalize the eigenvectors of h. (Of 

An v 

course the reciprocal property for the eigenstates of h also holds.) 

We remark that, depending on the specific form of M(x 2 ), one of the 
zero modes may not appear in the sum fllSD . In the Callan and Harvey 
mechanism, the defect is step-like, what gives rise to only one zero mode, 
with exponential localisation |3J . This issue will be made more explicit when 
considering concrete expressions for the mass function. 
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We then introduce the expansions ( |T5D into the Euclidean action (||), and 
see that it becomes a sum over actions corresponding to two dimensional 
fermions, labeled by the index of the eigenvalue 

S = sf + 4 0) + E / d 2 x^ n \x)(d + \ n )^ n \x) (18) 

where and S$ are the actions corresponding to the left and right zero 
modes, namely 

= Jd 2 x^\x)d^\x) , Sf = Jd 2 x^\x) d^\x). (19) 

As the fermionic integration measure decomposes into the infinite prod- 
uct: 

x Y[[D^ n) V^ in) } (20) 

it is evident that the fermionic determinant will be a product of two dimen- 
sional Euclidean determinants, one for each n. This yields an effective action 
which decomposes into an infinite sum: 

r = rf 3 + rg } + ]T (21) 

where 1^ and are the effective actions corresponding to the left and 
right zero modes, respectively: 

rf = -lndet(^) , rg J = -lndet(^), (22) 

and r*-™-* is the effective action for a massive two dimensional fermion with 
mass A„, 

rW = -lndet(j/+ A n ) . (23) 

Note that the sign of A n , which is chosen by convention, is irrelevant to 
the result of these two dimensional determinants. Namely, as the spacetime 
dimension for these modes is even, the mass term does not break parity and 
thus the fermionic determinant cannot depend on the mass sign. 

We now discuss, also for the case of a rectilinear defect, the relevance 
of the mass profile to the actual properties of the different contributions to 
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the effective action. Rather than considering defects of the most general 
possible kinds, we plan to study domain walls of a localized type, namely, 
defects having a typical width A. This width is defined by the extension of 
the region around a zero of the fermion mass, where it has an appreciable 
variation. The mass is assumed to have an approximately constant absolute 
value A outside the defect. This kind of behaviour is displayed in Figure 1. 
The analysis of the different contributions to the fermionic effective action 
has to take the finite width of the defect into account. To incorporate the 
localized modes, we expand the functional form for the mass around the 
defect in a Maclaurin series keeping a few terms, and then solve for the 
eigenmodes of h and h corresponding to that expansion. This procedure 
will only make sense if the modes so obtained are concentrated in a region 
contained in the band of width A. So we shall only add in the sum over 
modes fl2"T| ) those having a spatial dispersion smaller than A. This sum will 
give an approximate estimation of the contribution of the localized modes to 
the effective action. Non-localized modes correspond to Dirac fermions that 
live effectively in 2 + 1 dimensions, and shall have a smaller contribution to 
the effective action than the localized modes, as a simple argument shows: 
The action for a non-localized mode is given by the corresponding term in the 
The property that distinguishes it from the action for a localized 
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mode is that the former corresponds to an unbounded state of the operators 
h, h, while the latter is due to a bounded state. These operators can be 
written more explicitly as, 

h = -d 2 + M 2 (x 2 ) -d 2 M{x 2 ) 

h = -d 2 +M 2 {x 2 )+d 2 M{x 2 ) (24) 

which, far from the defect (for a defect like the one of Fig.l) may be replaced 
by 

h ~ h ~ -d 2 + A 2 , (25) 

since the derivative of the mass obviously tends to zero. We conclude that 
non-localized modes have an action, which has to be larger or equal than A, 
the absolute value of the mass outside the defect. 

If the mass function is approximately linear near to the center of the 
defect, we may keep only the first term in a Maclaurin expansion, arriving 
to the operators a and a) 

a = d 2 + M(x 2 ) = d 2 + M'(0)x 2 
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a t = -d 2 + M(x 2 ) 



-d 2 + M'(0) x 2 



(26) 



in terms of which we may define a and a\ which are identical to the usual 
(harmonic oscillator) creation and destruction operators: 



a = 



at = 



y/2M'(0) y/2M'(0) 



a, a 



1 . 



(27) 

We note that, as far as this linearization is reliable, we may of course write 
M'(0) in terms of A and A: 



M'(0) 



2A 

aT 



(28) 



The eigenvalues and eigenvectors of h and h are immediately found. We 
see that there is a zero mode for a but none for at , 



a</)o(x 2 ) = =>• O (^2) = A" e 2 



i|Af'(0)|xl 



(29) 



with Ao a normalisation constant. As the Hamiltonians fa and h are like 
number operators, the eigenvalues \ n = ^2|M'(0)|n are evenly spaced. 
Then the effective action becomes 



r = rf + r(n) 



n=l 



where 



p(n) 



Trln 



^+ ^2|M'(0)|n 



(30) 



(31) 



In order to determine the number A^ of modes that we may include in the 
sum above, we recall that for the harmonic oscillator states corresponding to 
the label n, the dispersion is 



a 



(n) 



rH-j 



(32) 



The condition that modes are localized inside the defect becomes = A, 
so 



N = [M'(0)A 2 



1 



(33) 
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As N > 0, we see that 



or 



M'(0)A 2 > 



A x A > 



(34) 



(35) 



which must be interpreted as a sufficient condition to fulfill in order to have 
at least one localized state. Its physical interpretation is quite transparent: 
the thinner the defect, the bigger should the mass jump be in order to localize 
a mode inside the defect. 

We can then exhibit the contribution to the effective action of these 
modes. The zero mode contribution is simply the effective action for a 
chiral fermion in 1 + 1 dimensions, a result we recall, for example, from 



reference |L1 



— I d 2 xA, 

87T 



a&tiv ~ (<W + it 



pa 



-02" 



(36) 



where a is a regularization parameter, which must be fixed to the value 1 for 
the regularization to be gauge invariant. 

Regarding the massive modes contributions, we apply the result presented 
in the Appendix for the massive determinant in the quadratic approximation, 
with a mass given by A n , to obtain the n mode contribution: 



rW(A) 



where 



and 



(2tt) 2 ^ 

fW(*) = fW(*)(V 



<hk AJ-k)V^{k)A v {k) 



h h 



A? 



f^(k) = -<h -2-f(i + 



4A 2 



7T 



k 2 



2 In 



(37) 

(38) 
(39) 



Then the effective action reads 

T(A) 



(]lk AJ-k)V v {k){5, v -^)A v {k) 



(2*) 



k' 2 



(40) 
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where 

? v (k) = £f (n) (*0 (41) 

n 

If a low momentum expansion for H n ) is used, then we have for T v the 
following expression: 



2 M'(0)A2-i 1 



n=l 



Moreover, for cases where the condition A A >> 1 (the number of local- 
ized modes is >> 1) is satisfied, the sum over n can be approximated by a 
logarithm, leading to 

r„ ~ f , x , [ln(AA) +7] A; 2 , (43) 

127r|M'(0)| L V ; n ' V ; 

where 7 is Euler's constant. 

In order to understand the nature of the error involved in the assumption 
of taking a linear profile for the mass near the defect, we recall the the 
expression for the error in a Maclaurin expansion to first order may be written 

as 

E(x 2 ) = \m"{9)xI (44) 

where < 6 < x 2 , and we assumed that M has a continuous second deriva- 
tive. A rough estimate of the value 5 of the radius of the interval in which 
we may use the linear approximation is by imposing the condition of having 
an error much smaller than the value of the linear term, which leads to 

5 « 2» , (45) 
\M"{9)\ 

where M" reaches its maximum at 9. 

We have considered so far only the case of one single defect. To see how 
the results for this case may be extended to the situation of having more than 
one defect, we note that our definition of localized modes assures that modes 
attached to different defects shall be, to a very good approximation, orthog- 
onal to each other. Then the effective action will be the sum of the effective 
actions corresponding to the modes in the different defects. This shows a 
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posteriori the utility of the truncation of the sum over modes according to 
their dispersion: it ensures the additivity of the corresponding contributions 
to the action. 



2.2 Defects of arbitrary shape 

We shall consider here the case of a wall having the shape of an arbitrary 
curve lying on the plane. For a curve C defined by a parametric equation 
(H), we introduce the (normalized) tangent, normal and binormal vectors, 
denoted by e 1? e 2 and e 3 , respectively To that end, it is better to use 
a natural parametrisation for the curve, defined in terms of the arc length 
parameter s 

s -> x(s) , s e [0, L] (47) 

where L is the total length of the curve. The tangent e\ is then obtained 
as a derivative of the parametric equations of the curve with respect to the 
length parameter 

ei = —S(s) . (48) 

It is natural to define an orientation for the curve: we shall assume that s 
increases when the region of positive mass is left to the right and the one 
of negative mass is to the left of the curve. The normal can, as usual, be 
defined by 

g 2 = ± \\ d ^( s ) u-i d ^( s ) (49) 
ds ds 

with a choice of sign which we fix by the convention that the normal should 
point from the side of positive to the one of negative mass. As the curve is 
plane, however, we can do more simply than that and define the normal as 
follows: 

ei = e**e{ (50) 

which is explicitly orthogonal to e\. It is immediate to realize that e*3 can be 
taken as a constant and pointing in the third direction, with a sense such that 



2 Of course, being the curve plane, it is torsionless (e 3 is constant), or, equivalcntly, the 
osculating plane coincides with the plane x% = for all the points on the curve. 
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dete* = +1. We now proceed to define coordinates in a neighborhood of 
the defect. Through each point on C, we draw a curve tangent to the normal 
vector, and another one tangent to the binormal. Of course the condition 
of being tangent only defines the curves in a small neighborhood, but this 
is enough for our purposes. We then perform a Lie dragging of the curve C 
along those fields. This defines a complete set of coordinate lines around the 
curve of the defect. 

Denoting by u±, u<i and u 3 the parameters labeling points along each one 
of the integral lines, we may then write the action corresponding to the region 
around the defect as 

Sc = J d 3 u det[ej]* [7XA + M( Ul )] V (51) 

where the spin connection term vanishes because the coordinate system has 
been defined by dragging. We then note that det[e*] = +1, that e l a di is the 
directional derivative along the integral curve a, and that e l a Ai is the gauge 
field component along the same direction. In order to arrive to a situation 
entirely similar to the one of the case of the rectilinear defect, we need to 
chose a convenient gauge. The analog of the condition A 2 = for that case 
would be here to use a gauge such that the component of A along the normal 
to the curve vanishes. We note that a condition like that does not forbid the 
presence of a non vanishing magnetic flux through the surface enclosed by 
the curve, since that requires only the existence of a non vanishing tangent 
component for the gauge field. We then arrive to the expression for the action 

S c = J d 3 u * [7 1 £> Ul + j 2 d U2 + ^D UA + M(u 2 )] V . (52) 

We conclude this section with an example of an application: the case of 
a circular defect. We assume the defect is a circle of radius R, and that the 
mass is positive inside the circle. The convenient choice of coordinates for 
this case corresponds to using polar coordinates. 

With the natural parametrisation, in terms of the arc length s, we have: 

x(s) = R( cos(^),sin(^) ) , < s < 2irR . (53) 
R R 

The tangent and normal vectors are, of course, 

ei = (-sin(-^),cos(-^)) 

e 2 = (cos(|),sin(|)). (54) 
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Local coordinates are defined as follows: u\ = rd, where r is the radial (plane) 
distance to the origin and 9 the polar angle (of course, u\ reduces to s when 
r = R). u 2 is equal to r — R and u 3 = x 3 . In terms of r, 9 and x 3 , the action 
)2f) for a defect of width A becomes 



S c = [ R+A drr d9 [ dx 3 ^ \-j\de + ieA e ) + 7 2 <9 r + j 3 D X3 + M {r - R) 
Jr-a Jo J lr 

(55) 

We may now consider, for example, a gauge field which corresponds to a 
static uniform magnetic field B(x) through the surface. Near to the defect, 
we may take 

A > = 2^ < 56) 
where $ is the total flux $ = BttR 2 . Note that, as we are only considering 
the fields near to the defect, the dependence of A e on the radius may be 
ignored, evaluating it at R. This will hold true as long as the radius of the 
defect is sufficiently large. The action for the chiral zero mode localized on 
the defect will then be 

Sl = R J q 27T d9 J dx 3 ^f\^i\d e + ze^) + l3 d 3 )^ ] (57) 

which corresponds to chiral fermions in the presence of a constant gauge field. 
This problem is equivalent to the one of free fermions with twisted boundary 
conditions, and the dynamics of the zero mode is then periodic in the flux. 
Non zero modes might also be incorporated, but only if their dispersion is 
small enough as to ignore de dependence of the gauge field on the radial 
direction. 



3 Fermionic current density 

In this section we present an evaluation of the expectation value of the 
fermionic current in the presence of defects, and under the influence of an 
external electric field. From what we said in the previous section, it should 
become evident that what we need is just to study the case of one single 
defect, and then the more general situations are easily understood from the 
knowledge of the results that follow for this case. We shall consider the vac- 
uum expectation value ^(^1,^2,^3) (/•* = 1,3) of the fermionic current in 
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the presence of an external gauge field such that F 2fJ/ = 0. That expectation 
value may of course be written as 

j»(x) = -etr [<y^S(x,x)] (58) 

where S(x, y) is the fermionic propagator, the inverse of the operator defining 
the quadratic form in the fermionic action: 

S af 3(x,y) = {x,a\V~ l \y,l3) (59) 

where 

V = (a+ <f)V L + (a f + $)V R . (60) 
The simple algebraic identity 

V- 1 = (Dtp)-ipt (ex) 

leads to 

V- 1 = (h- f)- x V L {ot- <f) + (h- f)- x V R {a- f) . (62) 

The x 2 dependence of the propagator is easily extracted by using the com- 
pleteness of the eigenfunctions of h and h, obtaining 

oo 

S a f3(x,y) = 4>n{x 2 )4>\ n {y 2 ){x 1 , x 3 , a\(\ 2 n - f) ~ l \ n V L \y x , y 3 , (3) 5 nm 

n=0,m=l 

oo 

- ( f ) n(x2)(f)l n (y2)(xi,x 3 ,a\(X 2 n - fy 1 fiV R \y 1 ,y 3 ,P)5 nm 

n=0,m=0 
oo 

+ ^n{x 2 )(f) ] m {y 2 ){x u x 3 ,a\{\ 2 n - ^^KPrIvi, y 3 , P)S nm 

n=l,m=0 
oo 

- (t>n{x 2 )$ m {y2) (Xl, X 3 , a\ (A*- f)' X $Pl\V\, VS, P)^nm ■ (63) 

ra=l,m=0 

Whence, the expectation value of the current becomes 

3im{x) = -e(/) 2 (x 2 )tr[>y ll (x 1 ,x 3 \ ^VrIx^ x 3 )] 

oo 

+ e J2 {4>l{x 2 )tr{-i^{xi,x 3 \(\ 2 n - f)' 1 $Vr\x 1 ,x 3 )] 

n=l 
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+4>l(x 2 )tr[j tJ ,(x 1 , x 3 \(Xl- 1 $V L \xi,x 3 )} 

-(pn(x 2 )^n(x2)tr['J^(x 1 ,X 3 \(\ 2 n - f) _1 \ n V L \ X 1 , X 3 )) 

-0 n (x 2 )0„(x2)tr[7 M (xi,x 3 |(A^- 2 )~ 1 \ n V R \x 1 ,x 3 }]} (64) 

The first term on the rhs is the contribution to the current coming from the 
chiral mode. Denoting by j„(x) and j^f 1 the vacuum expectation values of 
the two dimensional chiral and vector like currents, respectively, we see that 
the above result for the current may be rewritten as 

j^(x 1 ,x 2 ,x 3 ) = (f>l(x 2 )j^(x 1 ,x 3 ) 

1 oo 

- 9 E + #fcO)«W + - & 2 ))e M „] . (65) 



n- 



We will now apply the result (^) to the determination of the current due 
to a constant electric field E along the defects. It is evident from ( p5| ) that the 
x 2 dependence of the current is completely determined by the functions <p n 
and <j) n . For example, the contribution carried by the chiral zero mode shall be 
concentrated around the defects with a density 0q, which is a Gaussian for a 
linear defect. A step-like defect will in turn have an exponential localization. 
The ^-dependent factors for the n^O modes yield contributions which are 
more extended in x 2 , as the dispersion in x 2 grows with n. 

On the other hand, as the external electric field is assumed to be constant 
in time and uniform in space, there only remains to know the x\ dependence, 
which of course may appear in the expectation values of the two dimensional 
currents j£ and It can be easily checked that, within the linear response 
approximation, the expectation value of every massive current jj?' is zero, 
due to the existence of a gap | (the "mass" A n ). The contribution of the 
chiral zero mode is obtained from the parity- conserving part of its effective 
action, and yields exactly one half of the contribution corresponding to a 
massless Dirac fermion: 



tf( Xl ) = ^EL , (66) 



2 



where L is the length of the system in the direction x\. 



3 The conductivity of a massive fermion in 1 + 1 dimensions is proportional to the value 
of r(0), which vanishes for mass different from zero. It is however non-zero for a massless 
fermion. 
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We conclude that the three dimensional current is: 




where, as we have already mentioned, the form of <fil(x 2 ) depends on the 
particular shape of the defect. 

4 Conclusions 

Motivated by the fact that Dirac fermions with step like masses appear as 
effective descriptions of various systems in particular in Condensed Matter 
Physics, in this paper we have provided a formalism to study the general case 
in which the mass of the Dirac fermions changes sign along an arbitrary curve 
in the spatial plane, and with an arbitrary profile near the regions where the 
sign change is observed. 

We showed that the dynamics of the fermions can be studied in such 
a way that the contributions from the localized zero modes on the defects, 
and the ones of the higher no-zero modes can be clearly distinguished. We 
provided a criterion to evaluate which are the most relevant contributions 
depending on the profile of the mass near the region where it changes sign. 

We applied this formalism to the simplest case of a rectilinear domain wall 
and evaluate the vacuum fermionic current in the presence of an external field, 
to understand the consequences of having defects to the spatial distribution of 
the fermionic current. We studied also a circular defect with a magnetic field 
piercing the enclosed surface and saw that the effective theory for the zero 
mode is equivalent to the one of fermions with twisted boundary conditions 
in a circle. 
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Appendix: Two-dimensional massive fermionic 
determinant in the quadratic approximation 

We present here, for the sake of completeness, the calculation of the mas- 
sive fermionic determinant in two dimensions, in the presence of an external 



Abelian vector field [12j A^, in the quadratic approximation (namely, the 
first non-trivial term in the usual perturbation expansion). 

The generating functionals Z(A) and T(A), corresponding to massive 
Dirac fermions in the presence of an external Abelian vector field are 
defined by 

Z{A) = exp[-T(A)] = J VipVip exp{- J d 2 xijj(ft + ie A + m)i/j} 

= det($ + ie ft + m) . (68) 
Thus the effective action F is then given by 

T(A) = -lndet(^ + ze A + m) = -Tr ln(# + ie A + m) , (69) 

where the trace operates over both functional and Dirac space. In the 
quadratic approximation, T becomes 



r(A) 



--e 2 Ti 
2 



(70) 



+ m)- 1 4($ + m)- 1 4 
which, in momentum space representation, has the equivalent expression 

r(A) = \J ^- 2 A,(-k)f, u (k)A u (k) . (71) 

Here the tilde denotes the Fourier transform of a function of x, namely, 

Mk) = 



d 2 xe~ ik ' x A fl (x) 



d 2 x e~ 



ik-x 



(72) 



Fourier transforming in ([70]), we can give a more explicit expression for V as 
follows 



— e 



d 2 p 



tr 



1 



&+ $) + m i p 1 + m 



-7k 
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Gauge invariance implies the transversality of f (A;), which allows we to write 
r M „ in terms of only one independent scalar function T 

T(k) = f(fc 2 )(^-^). (74) 



We follow a procedure analogous to the one followed in [T(J for the mass- 
less case, to evaluate T using dimensional regularization. Thus, in d dimen- 
sions, after taking the trace on both sides of the previous tensorial equation, 



(d-l)f(A; 2 ) = 2e 2 (2-d) J 



d d p {k + p) ■ p 

(2Ti) d [(k + p) 2 + m 2 } (p 2 + m 2 ) 



-2 P 2 m 2 d ( ' AJL (7<\ 

J {27i) d [{k + p) 2 + m 2 }{p 2 + m 2 ) ' { ' 

Letting d — > 2, the first term on the right hand side of flT5|) contributes 
only through its divergent part, which is identical to the one for the case 
m = 

\^ {2 - d) J (ty(p 2 lrn 2 ) 2 = h' (?6) 

The second is regular at d = 2. Adding the contribution of both terms, we 
see that 

f (A; 2 ) = - - 4e 2 m 2 / -^j-, r= ^ . (77) 

V ; 7T J (2vr) 2 [(k + p) 2 + m 2 }(p 2 + m 2 ) V ; 

The integral over p is evaluated by using the standard Feynman parameter 
trick, to obtain 
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r , > e " I , r^™ 2 4m 2 . i _ 



4m' 
A; 2 



^+1 



i 



In a low momentum expansion, the leading contributions are 



2\2 



e^ijf_ (k 
it \ 6m 2 90m 4 



(78) 



(79) 
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Figure captions 

Figure 1 The typical mass profile considered, showing the definitions of 
the defect height 2A and width A. 
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